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I. INTRODUCTION 



Seeking out new physics beyond the standard 577(3) x 577(2) x U(l) model of particle 
interactions is the goal of the experimental program at the forthcoming Large Hadron Col- 
lider (LHC) at CERN. The discovery, and study, of such new physics will depend on our 
ability to calculate a wide variety of processes in the gauge theories making up the Standard 
Model. While computations of tree-level scattering amplitudes are a first step, the size and 
scale- variation of the strong coupling constant imply that a basic quantitative understanding 
must also include the one-loop amplitudes which enter into next-to-leading order corrections 
to cross sections. Such corrections are required to build a theoretical base for a program of 
precision measurements at hadron colliders [1]. Precision measurements at the SLAC Linear 
Collider (SLC) and CERN's Large Electron Positron (LEP) collider have proven the power 
of such a program in advancing our understanding of short-distance physics. 

Recently, there has been tremendous progress in understanding and computing gauge 
theory scattering amplitudes, stimulated by Witten's proposal [2] of a weak-weak coupling 
duality between M = 4 supersymmetric gauge theory and the topological open-string B 
model, generalizing Nair's earlier description [3] of the simplest gauge-theory amplitudes. 

While there are still open questions about the precise computation of tree-level ampli- 
tudes in the twistor string theory [4], and open questions about the structure of the theory 
at loop level [5], the proposal has already led to several important insights into the structure 
of gauge-theory amplitudes. In particular, it allowed Cachazo, Svrcek and Witten (CSW) [6] 
to formulate a new set of vertices and rules for computing tree-level amplitudes in either 
supersymmetric or non-supersymmetric gauge theories. The vertices are off-shell continua- 
tions of the maximally helicity-violating (MHV) amplitudes, which were written down long 
ago by Parke and Taylor [7] — in amplitude form by Berends and Giele [8] and by Mangano, 
Parke and Xu [9] — and proven using recurrence relations by Berends and Giele [10]. (For a 
discussion of the connections between different tree-level pictures, see refs. [11, 12].) A num- 
ber of papers have applied these techniques to explicit calculations [13, 14, 15]. Recently, 
several authors have taken the first steps in applying such ideas to mixed QCD-electroweak 
amplitudes as well [16]. At loop level, earlier apparent complications [17] in the twistor- 
space structure of amplitudes were clarified through an explicit linking of CSW vertices 
to the unitarity-based method [18], and via the elucidation [19] and computation of the 
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holomorphic anomaly [20, 21]. 

Motivated by early puzzles about the twistor-space structure of one-loop amplitudes, two 
groups pursued calculations of the next-to-MHV (NMHV) class of n-gluon amplitudes in 
the Af = 4 super- Yang-Mills theory [22, 23, 24, 25, 26], primarily using the unitarity-based 
method [27, 28, 29, 30]. A recent development, enhancing the power of the unitarity method, 
is the observation [24] that box integral coefficients can be obtained from generalized unitar- 
ity cuts [23, 30, 31] by solving the constraints that multiple intermediate states be on shell. 
(The notion of 'generalized unitarity', as applied to amplitudes for massive particles, can be 
traced back to ref. [32].) As a by-product of these calculations, ref. [23] revealed a simpler 
representation for the NMHV seven-gluon amplitude than had been known previously, from 
either a recurrence-based computation [33] or the CSW rules [6] . The NMHV computation 
for arbitrary values of n [25], like that for n = 7, made use of infrared consistency equa- 
tions, which express the constraint that a one-loop amplitude must have the correct infrared 
poles [34] in the dimensional regulator e. Again as a by-product, the solutions to the equa- 
tions revealed a new class of representations of tree amplitudes. A particularly nice one [35] 
gave the natural generalization of the simple representation to the eight-point amplitude, 
and suggested a new type of on-shell recurrence relation for tree amplitudes. Britto, Cac- 
hazo and Feng wrote down [36] the corresponding recurrence relation, which yields compact 
expressions for NMHV tree amplitudes [37] . The same authors and Witten (BCFW) gave a 
very simple and elegant proof [38] of the relation using special complex continuations of the 
external momenta. 

The proof, which we review in section III, is actually quite general, and applies to any 
rational function of the external spinors satisfying certain scaling and factorization prop- 
erties. The generality of the proof motivated us to find a loop version of the recurrence 
relation. In this paper, we shall apply it to the finite one-loop amplitudes in QCD, which 
are purely rational. We will also argue that it can be applied to compute the rational terms 
in other, cut-containing, QCD amplitudes. We shall see that the application to general loop 
amplitudes in a non-supersymmetric theory like QCD requires that we address new issues, 
such as double poles in the complex analytic continuation, which are related to properties 
of these amplitudes in the limit that two momenta become collinear. 

Recurrence relations can provide a technique, complementary to the unitarity-based 
method, for the complete calculation of one-loop QCD amplitudes. The unitarity method 
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applies most easily to terms in the amplitudes that have discontinuities. Computation of 
these terms requires only knowledge of tree amplitudes evaluated in four dimensions. The 
unitarity method can also be applied to computation of rational terms, because the latter 
effectively have cuts at higher order in e [29, 30]; but doing so requires knowledge of tree 
amplitudes with two legs in D = 4 — 2e dimensions, which are harder to compute. An- 
other approach is to build up an ansatz for the rational terms, based on their collinear and 
multi-particle factorization properties. This approach has been successful for amplitudes 
with sufficiently simple structure [39, 40]; and in more complex cases, with assistance from 
other, numerical representations of the amplitudes [31]. However, it involves some guesswork 
and so it is not easy to systematize. An alternative, recurrence-based method, along the 
lines of the present paper, would represent a systematic method for constructing rational 
function terms in amplitudes directly from their poles. It thus holds forth the prospect of 
dramatically simplifying such calculations. 

This paper is organized as follows. In section II we define the notation used in this 
paper. Section III reviews the BCFW proof. Then in section IV we construct recurrence 
relations for the all positive-helicity one-loop amplitudes. Section V presents such a relation 
for the one-loop amplitudes containing a single negative-helicity gluon, and also displays 
new compact forms for the cases of n = 6 and n = 7, obtained via the relation. Finally, in 
section VI we give our conclusions and outlook for the future. 



II. NOTATION 



We use the trace-based color decomposition [8, 9, 41, 42] of amplitudes. For tree- level 
amplitudes with n external gluons, this decomposition is, 

AT({hM,a l })= MT a ^---T a ^)AT(v(l hl ,---,n h ")), (1) 

where S n /Z n is the group of non-cyclic permutations on n symbols, and j hj denotes the j-th 
momentum and helicity hj. The T a are fundamental representation SU(iV c ) color matrices 
normalized so that Tr(T a T b ) = 5 ab . The color-ordered amplitude A^ ee is invariant under a 
cyclic permutation of its arguments. 

We describe the amplitudes using the spinor helicity formalism. In this formalism ampli- 
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tudes are expressed in terms of spinor inner-products, 

ao = (r\i + ) = «-(*,>+(*i) , b'fl = o' + io = «+(*>-(*i), (2) 

where u±(k) is a massless Weyl spinor with momentum k and plus or minus chirality [42, 43]. 
Our convention is that all legs are outgoing. The notation used here follows the standard 
QCD literature, with [ij] = sign(kfkj) (j i)* so that, 

(ij)[ji] = 2ki- kj = Sij . (3) 

(Note that the square bracket [ij] differs by an overall sign compared to the notation com- 
monly used in twistor-space studies [2].) 

We denote the sums of cyclicly-consecutive external momenta by 

K, = K + + + fc?_i + *?, (4) 

where all indices are mod n for an n-gluon amplitude. The invariant mass of this vector is 
Sj...j = Kf...j. Special cases include the two- and three-particle invariant masses, which are 
denoted by 

Sij = -K-ij = (ki + kj) = 2ki ■ kj, Sijk = [ki + kj + k^) • (5) 

In color-ordered amplitudes only invariants with cyclicly-consecutive arguments need appear, 
e.g. Sj ; j+i and s^+i^. We also write, for the sum of massless momenta belonging to a set 

A, 

K = £ K ■ ( 6 ) 

Spinor strings, such as 



and 



(r|(a + 6)|r> = M[aj] + (*&)[&j] , 
<i + | (a + 6)(c + d) \j~) = [ia] (a~\ (c + d) \j~) + [ib] (b~\ (c + d) \j~) , (8) 



will also make appearances. 
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For one-loop amplitudes, the color decomposition is similar to the tree-level case (1) [44]. 
When all internal particles transform in the adjoint representation of SU(iV c ), as is the case 
for Af = 4 supersymmetric Yang-Mills theory, we have 

|n/2j+l 

A 1 n - loop ({k l ,h z ,a t })= E Gr n . c (a) A n . c (a) , (9) 

C=l <T£Sn/ Sn;c 

where |_^J is the largest integer less than or equal to x. The leading color-structure factor 

Gr n; i(l) =N C Tr(T ai ---T a "), (10) 

is N c times the tree color factor. The subleading color structures are given by 

Gr n;c (l) = Tr(T ai ■■■T a <- 1 ) Tr(T flc • • • T a "). (11) 

S n is the set of all permutations of n objects, and S n - C is the subset leaving Gr n;c invariant. 

The one-loop subleading-color partial amplitudes are given by a sum over permutations 
of the leading-color ones [27]. Therefore we need to compute directly only the leading-color 
single-trace partial amplitudes A n -i(l hl , . . . ,n hn ). 



III. REVIEW OF THE PROOF 

The proof of the BCFW recurrence relation [38] starts by introducing a parameter- 
dependent shift of two of the external massless spinors, here labeled j and /, in an n-point 
process, 

Aj - ► Aj — zXi , 

A; ^X t + z\ j: (12) 
where z is a complex number. This induces a shift of the corresponding momenta, 

p? -pfW=pf-|0'"IV|r>, 

P f ^ P ?(z)=p? + z -(r\r\i-), (is) 

which preserves their masslessness, p'j(z) = = pf(z), as well as overall momentum conser- 
vation. This shift also implies, 

^p l (z)=p l + z(\r)(r\ + \j+)(i+\). (14) 
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Define 

A(z) = A(p u . . . ,p j (z),p j+1 , . . .,pi(z), . . . ,p n ), (15) 

which is an on-shell amplitude evaluated at a particular set of complex momenta. When 
A is a tree amplitude or finite one- loop amplitude, A(z) is a rational function of z. The 
physical amplitude is given by A(0). 
Now consider the following quantity, 



2m 

where the contour integral is taken around the circle at infinity. If A{z) — > as z — > oo, as 
in the tree-level cases considered by BCFW, then the integral vanishes, and we can evaluate 
this expression as follows, 

= A(0)+ V Res^, (17) 

Z ' 2 = Zo Z 

poles a 

where the sum is taken over the poles at the z a of A(z). We can then solve for A(0), 

A(0) = - ^ Rf^. (18) 

Z — Z a Z 

poles a 

If, on the other hand, A(z) — > as z — > oo, the contour integral is equal to C^, and 
this result is modified to, 

A(0)=C oo - Y Res^. (19) 

L — ' z=z a Z 
poles a 

As explained in ref. [38], so long as A(z) only has simple poles, each residue is given by 
factorizing the shifted amplitude on the appropriate pole in momentum invariants, so that 
at tree level, 

A(0) = Coo + A K Z = z rs)^A- R h (z = z rs ) , (20) 

r,s,h r "' s 

where h — ±1 labels the helicity of the intermediate state. There is generically a double 
sum, labeled by r, s, over momentum poles, with legs j and I always appearing on opposite 
sides of the pole. The squared momentum associated with that pole, Kf... s , is evaluated 
in the unshifted kinematics; whereas the on-shell amplitudes Al and Ar are evaluated in 
kinematics that have been shifted by eq. (12), with z = z rs . To extend the approach to one 
loop, the sum (20) must also be taken over the two ways of assigning the loop to the pair 
(A L ,A R ). 
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The poles in z arise from poles in momentum invariants; because of the structure of 
multiparticle factorization, in general only single poles arise. The same is almost true of 
collinear factorization at one loop; only lone powers of spinor products — and hence only 
single poles — appear in the splitting amplitudes in all helicity configurations except (+ + +) 

and ( ). In these identical-helicity cases, the tree splitting amplitude vanishes; however, 

the one-loop splitting amplitude has the form [a b] / (ab) 2 or its conjugate [27], and double 
poles do arise. In no case does anything worse than a double pole arise. As we shall see in 
section V, the double poles alter the form of the recurrence relation. 

At tree level, it is always possible to find choices of reference momenta j and I such that 
there is no boundary term; and only single poles arise. For the finite one-loop amplitudes 
we will consider explicitly in this paper, life is not quite so simple, and we must confront 
each of these problems in turn. 



IV. ALL-POSITIVE HELICITY ONE-LOOP AMPLITUDES 



Let us begin with the all-plus amplitude, for which Chalmers and the authors [39, 40] 
wrote a collinear-based conjecture [54], 



where 



and 



H n = 



96tt 2 (12) (2 3) •••((n- l)n) (nl)' 
^] Tr_ ^ 

l<U<i2<«3<*4<n 

Tr_ 1^ = - Tr[(l — 75)^^^^^^^] 
= (M2) M3] (hk) [Mi] • 



(21) 



(22) 



(23) 



The conjecture was later proven by Mahlon [45]. The factor N p counts the number of states 
circulating in the loop, +1 or —1 for each adjoint representation bosonic or fermionic state, 
respectively. For Dirac quarks in the fundamental representation, the contribution to N p 
is —2/N c , so that for QCD N p = 2(1 — rif/N c ), where rif is the number of quark flavors. 
Because the corresponding tree-level amplitude vanishes for this helicity configuration, the 
one-loop amplitude (21) is finite, and is purely a rational function of the spinor products. 
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With a simple shift as given in eq. (12), it is easy to see that there must necessarily be 
a boundary term. Without loss of generality, we may take / = n. There are two cases to 
consider: when j and n are adjacent, and when they are not. Take the second case. Just as 
in the tree-level MHV case, two denominator factors are shifted, 

((ra-l)ra) -> {(n-l)n)+z{(n-l)j) , (n 1) -> {n 1) + z {j 1) . (24) 

This gives rise to two distinct single poles in z, a result that could have been obtained 
without knowledge of the answer, because the form of the denominator is fixed by general 
collinear- and multiparticle-factorization arguments. 

Most terms in H n are either unshifted, or shifted by only a single power of z. These 
terms will fade away as z — > oo. The only dangerous terms are those involving both $j and 
f n . However, 

— fcjfin = = ~~^jkrikji (25) 

so the terms proportional to z 2 are those with i 2 = j and 14 = n, 

H n = z 2 \n~) {r\t {j+l) ... {n - l) \n~) 

= -z 2 (j-\$ 1 .. <j _ 1) \n-) 2 (26) 

= -z 2 {r\$ (j+1) ... (n _ 1) \ n -) 2 . 

Both the numerator and denominator scale as z 2 , giving rise to a boundary term. In the 
other case, when j and n are adjacent (that is, j — 1 or n — 1), only one pole (one power 
of z) is generated from the denominator; but correspondingly, the numerator also scales as 
z, so that again we obtain a boundary term (whose form is somewhat more complicated, 
though). 

It is not clear how to derive the boundary term without knowing the complete answer, so 
any recurrence relation we obtain with one is probably not useful for deriving the amplitude. 
We could, however, still use it for proving a conjecture. 

Because tree-level amplitudes containing less than two negative-helicity gluons vanish, 
we know that the one-loop identical-helicity amplitude has no multi-particle poles, and that 
the only collinear singularities come from the tree-level (— ; ++) splitting amplitude. This 
in turn implies that A(z) in this case has only single poles, and accordingly we can write a 
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recurrence relation for the all-plus amplitude, 

j n—l 1 

A£\l\...,n + ) = Oo+ EEE^,-,^-.^-^,)^- (27) 

r=l s=j v=0 r-"S 

xAt:l(K- h s , (s + 1) + , . . . ,n+, . . . , (r - 1)+) 

(where j ^ l,n — 1). For convenience we define A$ = A^ ee and = A^° op . The hatted 
momenta are evaluated using the shift (12), with 

K r-s 

% = ^rs = ~ — 77"" — ~~j — ZTT • (28) 



n 



Using the observation that tree amplitudes with a lone negative helicity continue to vanish 
generically for n > 3, even when evaluated at complex momenta [24], only four terms survive 
in the sum, 

V + 1) + , -K-j^A^ii^, (J + 2) + , (./ 1) + ) 

+W— 4 0) ((j - i) + >i + > (j + 1) + , u- 2) + ) 

(" 2)\K+_ hn )Af(-K^ n , (n - 1)+ 

+ -^4 1 2i(2 + , . . . , J + , . . . , (n - 1) + , ^i)4 0) (-^i^ + , 1 + )- (29) 

Here the three-point amplitude, before permuting its arguments and shifting its momenta, 
is given by 

^■ 2+ ^—Mry (30) 

With the choice of shift we have made, the three-point amplitudes A^(j + , (j + l) + , — K~ j+1 ) 
and A^((j — 1) + , j + , —KJ^J) vanish [36], so that the final relation takes the form, 

n +) iN P (j-\ft u+1) ... {n - 1) \n-) 2 



96tt 2 (1 2) (2 3) • • • ((n - 2) (n - 1)) ((n - 1) j) (j 1) 
+ ^^^i(l + , (n- 2) + ) ti,n)4° ) (-ti, n) (n - l) + ,n + ) 

+ -^42l(2 + , ••.../• (n 1) + , *i)4 0) (-*n,l> * + > 1 + )' (31) 

where z = — ((n — 1) n) / ((n — 1) j) and z = — (In) / (lj) in the last two lines, respectively. 
We have checked this relation numerically through n = 15, and find that it is indeed satisfied. 
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Xj — > Aj — zXi — z n , v A n , 



Can we find a recurrence which eliminates the boundary term? This is indeed possible, 
so long as we choose a more general shift of the original spinors. In particular, consider the 
following shift of the spinors (j < I < n), 

A, ^ X t + zXj, (32) 

A n T" Z . . Aj, 

which also keeps them massless and conserves overall four-momentum. The ratio (nj) / (I j) 
allows z to be assigned a definite weight under spinor phase rotations. So long as we choose 
/ not adjacent to n, and j adjacent to at most one of / and n, then we find that A(z) — > 
as z — > oo, and we obtain a recurrence free of surface terms. (We must accordingly take 
n > 5.) The price is that more terms must be included. For 1 < j < I — 1 and I < n — 1 
(which requires n > 6), we obtain, 

7r^-42i((i + 1)+, I' / (/ 2)\K+ hl )Af\-K r _ 1;l , (I - !)+,/+) 

+ 2) + , r ., (/ ^; +1 )4 0) (-^m + i, * + , (' + 1) + ) 



K 2 



/• (" " 2) + , X n + _l,n)4° j (-^n-l,n, (n " 1) + , ™ + ) 



42i(2 + , .-.../• /• in- l) + ,K+i) A f\-K^ + , 1+ ) • ( 33 ) 



K 2 ' 

The respective values of the shift variable z in the four terms in eq. (33) are 

((1-1)1) ((I + 1)0 (lj)((n-l)n) (lj)(ln) 



(34) 



((/-l)j)' ((/ + l)j)' {nj)((n-l)j)' (nj)(lj)- 

The recurrence relation (33) is depicted diagrammatically in fig. 1. For n — 5, there is 
an analogous relation containing only three terms. We have also checked the relation (33) 
numerically through n = 15. 



V. ONE-LOOP AMPLITUDES WITH A SINGLE NEGATIVE HELICITY 

Let us consider next n-gluon QCD amplitudes with a single negative-helicity gluon, and 
the remainder of positive helicity. These amplitudes were computed by Mahlon for an ar- 
bitrary number of external legs [45], using a loop- level generalization of the Berends-Giele 
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1-2 




n-2 



FIG. 1: Terms in the recursive expression (33) for An\l + ,2 + , . . . ,n + ). The vertices labeled by a 
T are trees, and the ones labeled by an L are loops. 

recurrence relations [10]. In this section, we wish to reproduce these amplitudes, but in 
a simpler form, following instead the BCFW recursive methodology. As we shall see, the 
essential difference between the tree amplitudes and the one-negative-helicity loop ampli- 
tudes is the appearance of double poles. As a result of their appearance, a straightforward 
application of the BCFW recurrence misses single pole terms underlying the double poles. 
We will find a modification which can be written, for the case at hand, in terms of the soft 
functions which control the universal behavior of amplitudes as momenta vanish. In general, 
boundary terms may arise, but they can be avoided by an appropriate choice of a pair of 
legs (j, /) to shift. 

It is instructive to first inspect the analytic properties of the five-point amplitude, as its 
features will serve as a guide to finding an all-n recurrence relation. The one-loop five-point 
amplitude with a single negative helicity leg was first calculated using string-based methods 
and is given by [39, 46], 



A?>(1- 2+ 3+4+ 5" 



96tt 2 (3 4> 2 
Now consider the effect of the shift, 



[2 5f 



+ 



(14) 3 [4 5](3 5) (13) 3 [3 2] (4 2) 



[12] [51] (12) (2 3) (4 5) 2 (1 5) (5 4) (3 2f 



(35) 



Ai 
Ai 
A 2 
A 2 



Ai, 

Ai — zA 2 , 
A 2 + z\i , 
A 2 . 



(36) 
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Applying it to eq. (35) gives the shifted amplitude, 



4^) 



N„ 



1 



96tt 2 (3 4) z 



[2 5] E 



+ 



;i4) 3 [4 5] (3 5) 



[12] ([51] -z[5 2]) 1 (12) ((2 3) + ,(13)) (4 5) 2 
(13) 3 [3 2] ((4 2) + ,(41)) 



(37) 



(15) (5 4)«3 2) + ,<31)) 2 . 

Since Af\z) vanishes as z — > oo, the choice of shift (36) does not generate a boundary 
term at z — > oo. Other shifts are also possible, but not all are as useful; some do generate 



boundary terms. Our shifted amplitude (37) has a new feature: the appearance of a double 
pole in z. The double pole implies that the BCFW recurrence relation cannot be applied 
directly, because it relied on the amplitude having only single poles, as is the case at tree 
level. Nevertheless, we will use eq. (37) to track down the required modifications. 
To proceed with our investigation, rewrite eq. (37) as a sum over pole terms, 



96tt 2 (3 4) z 



[2 5f 



;i4) 3 [4 5] (3 5) 



[12] ([51] -,[5 2]) (12)((23) + ,(13))(45) 2 
_ (13) 3 [3 2] (4 3) (21) 

" (15) (54) (31) ((32) + , (3 l)) 2 
(13) 3 [3 2] (41) 



(38) 



(15) (54) (31) ((32) + ,(31)). 

where we have used the Schouten identity to expose the single pole sitting under the double 
pole. Taking z — 0, we obtain a form of the amplitude whose terms we wish to map onto 
terms of a recurrence relation, 



A«(z = 0) = i 



96tt 2 



[2 5f 



(3 AY [12] [51] 



+ 



<14> 3 [4 5] (3 5) 
(3 4) 2 (12) (2 3) (4 5) 2 



+ 



(13) 2 [3 2] (41) 



(13) 2 [3 2](21) 
(3 4) (1 5) (4 5) (3 2) 2 (3 4) 2 (1 5) (5 4) (3 2) 



(39) 



The appearance of a double pole in eq. (38) may seem puzzling at first sight. It can be 
understood from the structure of the one-loop three-vertex with identical helicities, used for 
obtaining one-loop splitting amplitudes [27]. When the momenta of legs 1 and 2 become 
collinear, the vertex becomes [47, 48] 



4 1) (l+,2+,3+) = -i 



N n e 



(+) 

3 



(fa - k 2 



96tt 2 



V2 



.(+) J+) 



- 2 



fa-e^fa-e^- 
K\ 2 



—% 



. N p [12] [23] [31] 



96tt 2 



K 2 



(40) 
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FIG. 2: Terms in the recursive expression for A^\l , 2 + , 3 + , 4 + , 5 + ). 
It is useful to expose the kinematic pole, so we define 

2+ 3 + ) = -r^ 3 {1) (l + , 2\ 3+) . (41) 

K 12 

The propagator of the factorized leg provides a second factor of l/Kf 2 , leading to the 
observed double poles in the amplitude. These poles make their explicit appearance in 
eq. (35) as spinor products (2 3), (3 4) and (4 5) appearing quadratically in the denominator. 
(They can also be understood as arising from those collinear singularities, with k a parallel 
to k b , and k a ~ x{k a + kf,), which involve the one-loop splitting amplitude denoted by [27] 
Split^ loop (x,a + ,6+) = -tN p /(9Qn 2 ) x y/x(l - x) [ab] / (ab) 2 .) 

We turn now to the recursive construction of the five-point amplitude. Our input expres- 
sions are the one-loop vertex (41) as well as the one-loop four-point amplitudes 



4 1, (l+,2+,3+,4+) = -AMIl, (42) 
4 v ' ' ' ; 96tt 2 (12) (3 4) ' v ' 

Afd-XXA-) = .Ar#P^. (43) 

4 v ' ' ' ; 96tt 2 [12] (23) (34) [41] K J 

Although the symmetry is not manifest in the form we have written it, A^\l + , 2 + , 3 + , 4 + ) is 
invariant under cyclic permutations of its arguments. We shall also need the MHV n-gluon 
tree amplitudes [7, 9, 10], 

A$\l + , . . . , (j - 1) + , j~, (j + 1) + , . . . , (k - 1) + , k~,(k + 1) + , ■■■,n + )=i 

(44) 

and the three-point amplitude for (+- 1 — ), given in eq. (30). 
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The proof of the recurrence relations in ref. [38] leads us to expect a recurrence relation 
for the five-point one-loop amplitude of the form, 

^(l- 2 + ,3 + ,4+,5+) ~ 4° ) (i- ! -^i + 5 ) 5 + )^4 1) (2 + ,3 + ,4 + ,i, 1 - 5 ) 

n 15 

+ 4 0) (1-, 5+) * A?(2\ 3+, 4+, K+) 

21 15 



+ A? (i", ^ 2 + 3, 4 + , 5+) -L- 4 0) (2 + , 3+, -X 2 - 3; 



L 23 



+ 4 0) (i", ^23, 4 + , 5+) — ^- vf } (2+, 3 + , , (45) 

K JX 23) 



corresponding to the four diagrams illustrated in fig. 2. Following the discussion of sec- 
tion III, it is straightforward to evaluate each term. As we shall see, however, these terms 
do not yield the complete answer, because they miss a single pole lying underneath the 
double pole of the last term. However, an inspection of the discrepancy will allow us to 
determine an n-point recurrence relation. 

The first diagram, shown in fig. 2(a), corresponding to the first term in eq. (45), simply 
vanishes, 

D (a) = Af\l-, -K+, 5 + ) 4 1} (2 + , 3 + , 4+, = , (46) 

n 15 

because [36] 

4 0) (i-,-X+,5 + ) oc [MTi 5 )5] 3 oc (1-| (1 + 5) |5-> 3 = 0. (47) 

We may evaluate the next diagram, shown in fig. 2(b), corresponding to the second term 
in eq. (45), to obtain, 

£><»> = 4^i',-^r5 ! 5 + )^4 1) (2 + ,3 + ,4 + ,^ 5 ) 

-"■15 

_ N p (1 (-K 15 )f 1 [3 4] [2K 15 ] 



N P (1"| 


K 15 


\2-f 


96tt 2 wW(15) (5~| 


K 15 \ 


2-) 



X 7T- X 

-"■15 



[34](1-| 


K 15 | 


2-) 


(3 4) (2- 


I #15 


|2-> 



(48) 



where u = [K 15 2] and = (1K 15 ). Using 

(2-| K 15 |2"> = (2-| (1 + 5) |2"> + (1 _ |(1 ^ 5 5) |2 _ ) <1-| (1 + 5) |2-> = K 2 U , (49) 

and 



ww= (1~\K 15 \2~) = (15) [5 2] , (50) 
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diagram (b) of fig. 2 simplifies to 

D® = i Np 12 (51) 
96tt 2 (3 4) 2 [12][51] ' V ; 

This matches precisely the first term in our target expression (39). 

Similarly, it is not difficult to work out the diagram shown in fig. 2(c), with the result 

= 4 0) (i-,^ + 3,4 + ,5 + )-^4 0) (2 + ,3 + ,-X 2 - 3 ) 

-^23 



N p (14) 3 [4 5] (3 5) 
'96vr 2 (34) 2 (12) (23) (45)' 



^7^T2 7^T771^7^777V2 > \ bZ ) 



matching the second term of eq. (39). Finally, for diagram (d) we obtain the double-pole 
contribution very easily, 



N p (13) 2 [3 2](21) 



2 \2 
23 J 



= % i J X (53) 

96tt 2 (3 4)(1 5)(4 5)(3 2) 2 V ; 

matching the third term of the target expression (39). 

We now must confront the question of constructing the single-pole contribution in dia- 
gram (d) which underlies the double pole. This contribution must produce the fourth term 
in eq. (39). Notice that this term contains a factor of [3 2] / (3 2). Thus it should be dictated 
by the behavior of A^\l~ , 2 + , 3 + , 4 + , 5 + ) as the complex momenta k 2 and /c 3 become paral- 
lel. While for complex momenta the factor [3 2] / (3 2) can be singular as k<i • k% 0, for real 
momenta it is always nonsingular, because | (2 3) | = | [2 3] | = a/|2A; 2 • ks\. (Even for real 
momenta, the ratio does depend on an azimuthal rotation angle in this limit, in which k 2 and 
k 3 are rotated while holding k 2 + k s and k 2 ■ k 3 fixed.) As k 2 ■ k 3 — * 0, the factor [3 2] / (3 2) 
is subleading with respect to the leading collinear behavior captured by the one-loop split- 
ting amplitudes, Split+ loop (a;, 2+ 3+) oc [23] /(23) 2 and Split i" loop (a;, 2+ 3+) oc 1/(23). The 
usual discussions of universal collinear behavior of amplitudes limits [27, 47] do not extend 
to this level of accuracy. 

Instead we proceed by experimentation. To motivate a guess, we assume the single-pole 
term is related to the double-pole term (53) associated with K\ 3 = 0. Hence we try an 
ansatz which is equal to this expression, multiplied by some function which is suppressed as 
k 2 ■ k 3 — > 0. The additional multiplicative function should be dimensionless, and invariant 
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under phase rotations of spinors associated with all external states, and the intermediate 
state K.2%- We can use universal multiplicative "soft factors," which describe the insertion of 
a soft gluon s between two hard partons a and b in a color-ordered amplitude, to construct 
such a multiplicative function. The soft factors depend only on the helicity of the soft gluon 
and are given by [42], 

S {0 \a,s\b) = j^t, (54) 
(as) (so) 

S<°\a,8-,b) = -r^L (55) 
[a s\ [s b\ 

They are invariant under phase rotations of spinors associated with a and b, but not s. 
However, the product 

Kl 3 S^(a,s + ,b)S^(c,(-s)-,d) (56) 

is dimensionless, invariant under phase rotations of a, b, c, d and s, and suppressed as 
K% 3 — > (for suitable choices of a, b, c, d and s). 

Since s appears in two factors in eq. (56) carrying opposite helicity, it is natural to identify 
it with the on-shell intermediate momentum K 2 3- Choosing c = 3 and d — 2 produces the 
desired collinear behavior oc (2 3) for eq. (56). With a little more experimentation, we find 
that in this case, the color- adjacent legs 1 and 4 should be identified with a and b. Thus we 
arrive at the candidate single-pole term given by 



Df ] = Af\l-, Kn, 4 + , 5 + ) 5<°>(i, K+, 4)^ S(°>(3, -K„, 2) vf (2+, 3 + , -K+) . (57) 
Evaluating eq. (57) we have 



K 2 

21 23 



D?=i 3L 2 <13>-[32]<41> f (5g) 

1 96tt 2 (3 4) 2 (15)(5 4)(3 2) V ; 

so it reproduces the last term in the amplitude (39). 

We now take the value of the diagram in fig. 2(d) to be the sum of the single- and 
double-pole contributions, 



D (d) = D (d) + D {d) 



i- 



(13) 2 [3 2] (21) (13) 2 [3 2] (41) 



+ 



(3 4) (1 5) (4 5) (3 2) 2 (3 4) 2 (1 5) (5 4) (3 2) 



96tt 2 

which matches the last two terms in the amplitude (39). 



(59) 
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FIG. 3: The diagrams describing a recursive formula for A n ' (1~ , 2 + , . . . , n + ). Diagram (a) van- 
ishes, while diagram (d) needs to be summed over j. For j = 3 diagram (d) contains both single-pole 
and double-pole contributions. 

We have thus succeeded in finding a simple interpretation of the five-point amplitude 
(39) in terms of a recursive diagrammatic construction, 

4 1) (l",2 + ,3 + ,4+,5 + ) = + + + , (60) 

corresponding to the diagrams in fig. 2, after interpreting diagram (d) as a sum over single- 
pole and double-pole contributions. 

To generalize this construction to n legs, we assume that the term corresponding to 
the single pole underlying the double pole in K% 3 is given by the double-pole contribution, 
multiplied by the same function found to work for n — 5. The resulting recursion relation, 
shown diagrammatically in fig. 3, is 

4 1} (l-,2 + ,...,n + ) 

= A n l \{A\ 5+, . . . , n\ i" K+) ±- Af\2\ 3+, -K^) 

-^23 

+ E A n- j+ 2(U + U + 2) + , • • • , n + , i", K-J * Af{2\ 3+, . . . , j\ -K+J 

+ ^(4+, 5 + , . . . , n+, X 2 - 3 ) ? * \/ 3 (1) (2 + , 3+, -K+) 

x (l + iv 2 2 3 5( )(i,iv 2 3,4)5(°)(3,-^ 2 3,2)). (61) 

The hatted momenta are evaluated using the shift (36), where the value of z for the K 2 ...j 
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channel is, 



K 2 



2...j 



- d-|^|2->- ^ 
For self-consistency, amplitudes generated by this recurrence relation must vanish under 

the shift (36) in the limit z — > oo. Otherwise, there would be a boundary term not included 

in the recurrence relation. To show that the amplitudes generated by eq. (61) satisfy this 

criterion, we first observe that under the shift the explicit poles appearing in the amplitude 

will behave either as 1/z for the single poles or 1/z 2 for the 1/(K 2 3 ) 2 double pole. In fact, 

these are the only places z-dependence will arise. (Basically, the point of the recursive 

decomposition is to expose the z-dependence as pole factors multiplied by ^-independent 

residues.) To see this explicitly, first note that Ai is unaffected by the shift, while Ai is 

affected. In the jth term in eq. (61), consider the spinor product [la] for an arbitrary leg 

a (not equal to K 2 „,-). It looks as though it might develop z-dependence in the next shift. 

However, what actually appears is the hatted expression, 

K 2 



[ia] = [lBl Vi*,.>> [2fll 

( 2+ l ft(j+l)-lft(j+l)-n \ a 



(2 + \it {j+1) ..J 1 \ 


a') -(2+| t 




a-) 


(11 




2-) 



(63) 



(1"| W(j+iy.. n l 2 ~) 

Since Ai and A 2 are unaffected by the shift, the only possible source of z-dependence is from 
the inside Jjt^-j- But this term is proportional to [2 1], and thus unaffected (because the 
Ai shift is proportional to A 2 ). 

Similarly, the expressions on the other side of the K 2 ..^ pole which might appear to 
develop z-dependence are of the form 

<^< 2t >- ( r|£|2-) <16> 

(64) 

and so they also are easily seen to be unaffected by the subsequent z shift. An expression 
involving the intermediate hatted momentum, like {2K 2 ...j), can be re-written as 

(65) 

and then Jj£ 2 ...j can be split into a sum of momenta in the numerator, applying eq. (64) to 
each term. The denominator is also inert. Thus the only z-dependence in eq. (61) is from 
the explicit K\___- pole factors, and the expression vanishes as z — > oo. 
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(1-1 




-3 


6+) 


(11 




2-) 



(2"| 


$2- 3 


2~) 


(11 


$2... 3 


2-) 



We have checked the recurrence relation (61) numerically against the computation of 
Mahlon [45] through n = 15. The relation leads straightforwardly to compact expressions 
for higher-point amplitudes. For example, the six-point amplitude can be written as, 



A> ;1 (l-,2+,3 + ,4 + ,5+,6+) 



. N p 
'96tt 2 



-\3 



(11(2 + 3)16-)' (1-|(3 + 4)|2-) 



+ 



(12)(2 3)(4 5) 2 s 123 (3-|(l + 2)|6-) (3 4) 2 (5 6) (6 1) s 234 (5"| (3 + 4) |2~) 

[2 6f ( [2 3] [3 4] [4 5] [5 6] [3 5] \ 

[12] [61] S3 45 V(4 5) (5- 1 (3 + 4) |2") (3 4) (3"| (1 + 2) |6") + (3 4) (4 5) J 
(13) 3 [2 3](2 4) (15) 3 (4 6) [5 6] 



(2 3) 2 (3 4) 2 (4 5) (5 6) (61) (12) (2 3) (3 4) (4 5) 2 (5 6) 2 
(14) 3 (3 5)(1-|(2 + 3)|4-) 



(12) (2 3) (3 4) 2 (4 5) 2 (5 6) (61) 



(66) 



Similarly, in the seven-point case we obtain the form, 

A 7;1 (l- 2+ 3+ 4+ 5+ 6+ 7+) = if + if | flip , (67) 

where 

X(l, 2, 3, 4, 5, 6, 7) | flip = -X(l, 7, 6, 5, 4, 3, 2) , (68) 

and 



7 96tt 2 



-\3 



(11(2 + 3)|7-) 



(12) (2 3) (4 5) (5 6) s 123 s 456 (3"| (1 + 2) |7") 



g45 (l-|(2 + 3)|4-) (6+|(4 + 5)(5 + 6)|7-) 
[ (l-|(2 + 3)(4 + 5)|6+) + (4-|(5 + 6)|7-) 
(1-|(5 + 6)|7-) 3 



(12) (23) (34) (56) 2 s 567 (4"| (5 + 6) |7~) 
, [2 7] 3 (1 [3 6] + (6 + | (4 + 5) (5 + 6) |7~) 



[12] (3 4) (4 5) (5 6) [71] s 712 \2 L J (3"| (1 + 2) |7~) 

1 S45 (2+|(3 + 4)(5 + 6)|7-) 



2(3-|(l + 2)|7-) (6"|(7+1)|2- 
(13) 3 [23](24) (14) 3 (35)(1-|(2 + 3)|4-) 



(2 3) 2 (3 4) 2 (4 5) (5 6) (6 7) (7 1) (1 2) (2 3) (3 4) 2 (4 5) 2 (5 6) (6 7) (7 1) 

_1 (l-|(2 + 3)(4 + 5)|l+) 3 

2 (1 2) (2 3) (4 5) 2 (6 7) (7 1) (1" | (6 + 7) (4 + 5) |3+) (1"| (2 + 3)(4 + 5) |6+) 

(69) 

Although the amplitudes (66) and (69) are quite compact in comparison with the 
expressions appearing in ref. [45], they do contain spurious denominators of the form 
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(a - 1 (b + c) \d~) and (a~\ (b + c)(d + e) \f + )- Such denominators typically arise from reduc- 
ing pentagon integrals to lower-point integrals. It is of course no surprise that we generate 
such denominators, given the origin of the BCFW-type recurrences in the relations between 
coefficients of one-loop amplitudes. Similar spurious denominators also arise in the CSW 
construction of amplitudes from MHV vertices. In fact, it appears empirically that spurious 
denominators are necessary to obtain the most compact forms of amplitudes that contain 
many multi-particle factorization poles. 

VI. CONCLUSIONS AND OUTLOOK 

In this paper, we presented a recursive approach to obtaining rational terms appearing 
in non-supersymmetric one-loop amplitudes, based on the Britto, Cachazo, Feng and Wit- 
ten [36, 38] recursive approach to tree-level amplitudes. A striking feature of the proof of 
the method in ref. [38] is that it relies only on very general properties of field theory and 
complex function theory. It should therefore apply much more broadly than to massless 
gauge theory tree amplitudes. This motivated us to revisit the question of loop amplitudes. 

In this paper, as an initial application, we focused on the finite one-loop all-gluon QCD 
amplitudes which are purely rational: amplitudes with all gluons of positive helicity, and 
those with a lone negative-helicity gluon. These amplitudes had been obtained previously in 
refs. [39, 40, 45], and are therefore quite useful as a laboratory for studying the new features 
appearing at loop level. Two essential differences arising at loop level are boundary terms 
and double poles, both of which modify the form of the recurrence. The boundary terms 
interfere with obtaining a useful recurrence relation, but via an appropriate choice of shifts 
we can avoid these. From general field theory considerations, tree-level amplitudes cannot 
have double poles in kinematic invariants. However, at loop level double poles do arise, when 
using complex momenta. These double poles induce extra terms in the recurrence relations, 
which appear to be related to the universal functions describing soft gluon emission, at least 
for the case discussed here. 

In general, we envisage that the logarithmic parts of one-loop QCD amplitudes will be 
obtainable from unitarity or related methods, while the rational function parts will be calcu- 
lable using recursive methods. We can think of all-gluon QCD amplitudes as being composed 
of three distinct terms: an Af = 4 supersymmetric amplitude; a multiple of the contribution 
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of an M = 1 supersymmetric chiral multiplet; and a multiple of the contribution of a colored 
scalar circulating in the loop. (There is a similar type of decomposition for amplitudes with 
external fermions [50].) The first two terms can be computed from their four- dimensional 
cuts, for which in turn the computation of massless four-dimensional helicity amplitudes suf- 
fices. Many such amplitudes have already been computed, in both the M = 4 [22, 23, 25, 27] 
and M = 1 cases [28, 51, 52] 

The scalar term, which is non-supersymmetric, can also be computed via unitarity, but the 
computation would require tree amplitudes with the scalar propagating in the full D = 4 — 2e 
spacetime [29], which are more cumbersome for large numbers of external legs. We can, 
however, think of splitting the scalar contribution up further into those terms which can 
be computed from the four-dimensional cuts, and a remainder. Indeed, the cut-containing 
parts of the scalar loop amplitudes with two negative-helicity gluons and (n — 2) of positive 
helicity have already been computed [28, 53]. To the terms with discontinuities such as 
dilogarithms and logarithms, it is convenient to add certain rational terms which remove 
spurious (unphysical) singularities [46]. For example, the function ln[(— sx)/(— s 2 )]/ (s 1 — s 2 ) 2 
generically appears in one-loop amplitudes, and has a singularity as the momentum invariant 
si approaches s 2 , which is not associated with any physical factorization. This singularity 
can be removed by replacing 



where Li(r) = (Lo(r) + l)/(l — r),Lo(r) = m(r)/(l — r). Subtracting from the full amplitude 
all the cut- containing terms, and those rational functions linked to them by spurious sin- 
gularities, defines a 'pure-rational' remainder. This function may be computed by on-shell 
recurrent techniques as done at tree level by BCFW [38], and at one loop in the present 
paper. To analyze the factorization properties, that is the residues at the poles in z, it may 
prove useful to think of the pure-rational remainder as the difference between the scalar 
loop amplitude in D = 4 — 2e and that in D = 4. Terms in the loop integral which are 
proportional to the (— 2e) components of the loop momentum are well-behaved in the in- 
frared. For this reason, complications (discontinuities) in amplitude factorization which are 
associated with infrared divergences [47] should be absent from the appropriate differences 
of amplitudes. 

One may wonder whether one can push this approach further. Indeed, it seems quite likely 




(70) 
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that it will apply at tree level to a wide variety of field theories, including scalar field theories 
or gauge theories with massive fermions or scalars, and indeed even to spontaneously-broken 
gauge theories. Another question is to what extent it can be applied beyond tree level, 
and indeed beyond the rational terms at loop level as discussed above. At loop level, the 
amplitude will now have branch cuts in z; but in general, we may still expect it to fall off at 
infinity, allowing A(0) to be analyzed in terms of simpler (lower-point) amplitudes. We are 
optimistic that these techniques will open the door to many new computations of Standard 
Model scattering amplitudes. 
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